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1. Introduction 

Doppler-free spectra of the near- infrared N I transitions in the 2p 2 ( 3 P)3s 4 P — > 
2p 2 ( 3 P)3p 4 P° and 2p 2 ( i P)3s 4 P 2p 2 ( 3 P)3p A D° multiplets have been recorded 
by Jennerich et al. [1\ using saturated absorption spectroscopy, extending the pioneer 
work of Cangiano et al. [2] using a similar set-up. From the analysis of these spectra, 
Jennerich et al. [1] determined the hyperfine structure constants for the various J-values 
of the three multiplets involved, for both isotopes 14 N and 15 N. The isotope shifts in 
each multiplet have also been measured, revealing a significant J-dependence of the 
shifts. These authors recommended a theoretical investigation of the underlying cause 
of this unexpected phenomenon. Similar measurements, using Doppler-free saturated 
absorption spectroscopy, have been performed in atomic Fluorine [3], Chlorine [3] and 
Oxygen [5]. 

The present work presents a robust ab initio theoretical estimation of the relevant 
hyperfine structure parameters, using the ATSP2K package J6], based on the non 
relativistic multiconfiguration Hartree-Fock (MCHF) and configuration interaction (CI) 
methods. The minimum theoretical background is presented in section [2] for the 
hyperfine interaction. Section [3] is dedicated to the description of the electron correlation 
models. A refined calculation for the even-parity term is presented in section HJ 
A first comparison between theory and observation appears in section [5l revealing 
unexpected discrepancies. Relativistic effects are investigated in section [6] through the 
relativistic configuration interaction (RCI) approach, confirming serious problems in the 
experimental analysis as discussed in section [7J 

2. The hyperfine interaction 

The theory underlying MCHF calculations of hyperfine structure parameters can be 
found in [7J E]- Neglecting the relativistic effects, the diagonal and off-diagonal A 
and B hyperfine interaction constants are expressed in terms of the J-independent 
orbital (a/), spin-dipole (a s d), contact (a c ) and electric quadrupole (b q ) electronic 
hyperfine parameters defined by 
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and calculated for the magnetic component Ml = L and Ms = S. The first three 
parameters (PQ), ([2]) and ([3]), contribute to the magnetic dipole hyperfine interaction 
constant through 

Aj = AJ b + Af + A c j, (5) 

with 

d _ 1 ^ 3(L.5 ? )(L.J) - L(L + 1)(5 ? .J) 
Aj "2 G ^ s T Gsrf SX(2L-1)J(J+1) ' ( j 

Aj _ 6 M T 5J(J+1) ' (8) 
while the last one constitutes the electronic contribution to the electric quadrupole 
hyperfine interaction 

6(L.J) 2 - 3(L.J) - 2L(L + 1)J(J+1) 
B '-- G *Q b * L(2L-l)(J+l)(2J + 3) • (9) 

Expressing the electronic parameters a/, a s d and a c in atomic units (flQ 3 ) and [Li 
in nuclear magnetons (//jv), the magnetic dipole hyperfine structure constant A is 
calculated in units of frequency (MHz) by using = 95.41067. Similarly, the electric 
quadrupole hyperfine structure constant B is expressed in MHz when adopting atomic 
units (ctQ 3 ) for b q , barns for Q and G q = 234.96475. 

The electronic parameter governing the mass isotope shift of an atomic energy level 
is the S sms parameter [9] defined by 

TV 

S S ms = -(jLSM L M s \J2^i ■ Vj| jLSM l M s ) . (10) 

i<j 



3. Correlation models 

The multiconfiguration Hartree-Fock (MCHF) variational approach consists in 
optimizing the one-electron functions spanning a configuration space and the mixing 
coefficients of the interacting configuration state functions (CSF) [10] for describing a 
given term 

V(ryLSM L M s ) = J2^(^LSM L M S ). (11) 

i 

Efficient MCHF expansions are often built by allowing single and double excitations from 
a multireference set (SD-MR-MCHF). As far as hyperfine structures are concerned, suc- 
cessful applications of this method are found for light elements such as Li-like ions [9], 
Be I, B I, C II and C I PH G2], N I P^, O I [131 or Na I PS]. 
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For the even parity term 2p 2 3s 4 P, the configuration expansion is generated from 
single and double (SD) excitations from the 



multireference (MR) to increasing active sets (AS) of orbitals that are denoted by spec- 
ifying the number of orbitals for each /-symmetry. This multireference set captures the 
dominant correlation effects through a physical good "zero-order" wave function. We 
include only configuration state functions ($(7jLS 'M lM '$)) that interact with the mul- 
tireference, adopting the reversed orbital order, i.e. coupling sequentially the subshells 
by decreasing n and /. This technique reduces substancially the size of the MCHF ex- 
pansions while keeping the dominant correlation contributions [121 [16]. With the largest 
set of orbitals (10s9p8d6f3g) optimized through these calculations, the effect of higher 
excitations is investigated through configuration interaction (CI) calculations using the 
configuration state function set obtained by adding to the original SD-MR CSF list, the 
triple and quadruple excitations from the same multireference to smaller orbital active 
sets (up to 6s5p4d3f). For generating these lists, some limitations have been introduced 
by imposing the restriction that there should be at least 5 orbitals with n < 3 in the 
CSFs produced. The merging of the original SD-MR and the TQ-MR CSF lists is noted 
by the union (U) symbol. 

For the odd parity states 2p 2 3p A P° and 2p 2 3p 4 D° ) a similar strategy is adopted, 
using the following multireference set 



The total energies (E), the specific mass shift (S sms ) and the hyperfine interactions 
parameters {ai, a s d, a c , b q } defined in section [2} are reported in Tables [U [2] and EJ 
respectively. These tables illustrate the smooth convergence of the various parameters 
while improving the correlation model by increasing the size of the orbital active set. 
For the even parity term 2p 2 3s 4 P (Tabled]), triple and quadruple excitations from the 
multireference affect the ai and a s d parameters at the level of 2% while the contact term 
a c is much more sensitive (33%). For the odd parity 2p 2 3p A P° and 2p 2 3p A D° terms 
(Tables |2] and [3]), all hyperfine parameters reach a high degree of convergence, except 
the contact parameter that is strongly affected (around 30%) by TQ excitations. For 
the three terms, the last layer added at the SDTQ-CI level of approximation brings 
a negligible increment, except for the contact contributions (3%). The convergence 
patterns of the contact contributions, as well as the sensitivity for TQ excitations, are 
similar to the one found for the 2p 3 3s 5 S° term in O I |13j . 

4. On a larger Multireference set for 2p 2 ( 3 P)3s 4 P 

We investigate the reliability of the theoretical parameters by extending the 
multireference set. Amongst the three terms considered, we focus on the even parity 
2p 2 ( 3 P)3s A P one, realizing that it is one for which the hyperfine interaction parameters 



{2s 2 2p 2 3s, 2p 4 3s, 2s2p 4 } 



(12) 



{2s 2 2p 2 3p, 2p 4 3p, 2s2p 2 3p3d} . 



(13) 
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Table 1. SD-MR-MCHF and CI total energies, specific mass shift and hyperfine 
interaction parameters (all in atomic units) as a function of the orbital active set, for 
2p 2 3s 4 P. NCSF is the number of CSFs in the wavefunction expansion. 



AS 


NCSF 


E 




Ssms 


ai 


dsd 


a c 


\ 








SD 


-MR-MCHF 










HF 


1 


-54.032 


303 


-2.5647 


3.7993 


0.7599 


3.6055 


1.5197 


3s2pld 


149 


-54.116 


873 


-2.3529 


3.5643 


0.7310 


8.4493 


1.2085 


4s3p3dlf 


652 


-54.165 


015 


-1.3937 


3.7020 


0.7999 


9.2913 


1.2704 


5s4p3d2flg 


1626 


-54.189 


786 


-1.4252 


3.6990 


0.7837 


7.3833 


1.2072 


6s5p4d3f2g 


3082 


-54.196 


747 


-1.4254 


3.6862 


0.7713 


8.7173 


1.2108 


7s6p5d4f3g 


5020 


-54.199 


972 


-1.4238 


3.6869 


0.7723 


8.1898 


1.2281 


8s7p6d5f3g 


7113 


-54.201 


188 


-1.4238 


3.6879 


0.7751 


8.0753 


1.2248 


9s8p7d6f3g 


9572 


-54.201 


746 


-1.4237 


3.6878 


0.7752 


8.0941 


1.2208 


10s9p8d6f3g 


11728 


-54.201 


953 


-1.4234 


3.6884 


0.7743 


8.1090 


1.2231 


SDTQ-CI 


U 4s3p2dlf 


40685 


-54.203 


791 


-1.4282 


3.6412 


0.7668 


9.7166 


1.1502 


U 5s4p3d2f 


106472 


-54.204 


677 


-1.4427 


3.6170 


0.7624 


10.511 


1.1156 


U 6s5p4d3f 


210533 


-54.204 


967 


-1.4470 


3.6090 


0.7609 


10.780 


1.1049 



Table 2. SD-MR-MCHF and CI total energies, specific mass shift and hyperfine 
interaction parameters (all in atomic units) as a function of the orbital active set, for 
2p 2 3p A P° . NCSF is the number of CSFs in the wavefunction expansion. 



AS 


NCSF 


E 


$ sms 


ai 


a s d 


a c 


\ 


SD-MR-MCHF 


HF 


1 


-53.984 055 


-2.5786 


1.9609 


-0.3749 


0.0000 


-0.7498 


3s2pld 


583 


-54.062 052 


-2.2945 


1.9466 


-0.3787 


5.1320 


-0.7123 


4s3p2dlf 


3879 


-54.112 044 


-1.2732 


1.9752 


-0.4086 


0.5548 


-0.6912 


5s4p3d2flg 


10078 


-54.134 423 


-1.3540 


1.9653 


-0.3884 


0.3990 


-0.6708 


6s5p4d3f2g 


19200 


-54.143 169 


-1.3318 


1.9673 


-0.3839 


1.1424 


-0.6841 


7s6p5d4f3g 


31245 


-54.147 455 


-1.3170 


1.9709 


-0.3835 


0.7422 


-0.6879 


8s7p6d5f3g 


44096 


-54.148 983 


-1.3187 


1.9705 


-0.3837 


0.6826 


-0.6887 


9s8p7d6f3g 


59110 


-54.149 672 


-1.3185 


1.9702 


-0.3850 


0.7044 


-0.6864 


10s9p8d6f3g 


72070 


-54.149 921 


-1.3188 


1.9704 


-0.3843 


0.7449 


-0.6849 


SDTQ-CI 


U 4s3p2dlf 


124029 


-54.150 179 


-1.3136 


1.9710 


-0.3844 


0.8878 


-0.6848 


U 5s4p3d2f 


252690 


-54.150 330 


-1.3139 


1.9714 


-0.3844 


0.9168 


-0.6845 


U 6s5p4d3f 


459494 


-54.150 406 


-1.3138 


1.9717 


-0.3843 


0.9429 


-0.6845 



Aj are the most sensitive to TQ excitations due to the large contribution of the contact 
term. The following multireference (MR') 

{2s 2 2p 2 3s, 2p 4 3s, 2s2p 4 , 2s2p 3 3p, 2s2p 2 3s3d} (14) 

is selected, after a detailed analysis of the eigenvector weights obtained with the first 
approach. The results are reported in Table HI The comparison of the last lines 
of Table [1] (MR) and Table H] (MR') illustrates the global stability of the hyperfine 
parameters, the largest variation (3.6%) being observed for the a c contact contribution. 
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Table 3. SD-MR-MCHF and CI total energies, specific mass shift and hyperfine 
interaction parameters (all in atomic units) as a function of the orbital active set, for 
2p 2 3p A D° . NCSF is the number of CSFs in the wavefunction expansion. 



AS 


NCSF 


E 




Ssms 


ai 


dsd 


a c 


\ 








SD 


-MR-MCHF 










HF 


1 


-53.987 


051 


-2.5789 


3.9263 


0.7462 


0.0000 


1.4926 


3s2pld 


720 


-54.064 


819 


-2.2937 


3.8979 


0.7539 


5.1069 


1.4194 


4s3p2dlf 


4747 


-54.114 


835 


-1.2729 


3.9574 


0.8137 


0.5229 


1.3775 


5s4p3d2flg 


12197 


-54.137 


237 


-1.3541 


3.9575 


0.7769 


0.3613 


1.3421 


6s5p4d3f2g 


23072 


-54.146 


073 


-1.3259 


3.9674 


0.7698 


1.1034 


1.3472 


7s6p5d4f3g 


37372 


-54.150 


182 


-1.3163 


3.9688 


0.7661 


0.6698 


1.3768 


8s7p6d5f3g 


52558 


-54.151 


947 


-1.3186 


3.9688 


0.7658 


0.6258 


1.3800 


9s8p7d6f3g 


70266 


-54.152 


646 


-1.3188 


3.9678 


0.7692 


0.6472 


1.3730 


10s9p8d6f3g 


85436 


-54.152 


903 


-1.3191 


3.9682 


0.7688 


0.6864 


1.3686 


SDTQ-CI 


U 4s3p2dlf 


156938 


-54.153 


165 


-1.3139 


3.9695 


0.7691 


0.8307 


1.3684 


U 5s4p3d2f 


334710 


-54.153 


318 


-1.3142 


3.9706 


0.7690 


0.8600 


1.3678 


U 6s5p4d3f 


620614 


-54.153 


406 


-1.3138 


3.9711 


0.7689 


0.8848 


1.3677 



Table 4. SD-MR'-MCHF and CI total energies, specific mass shift and hyperfine 
interaction parameters (all in atomic units) as a function of the orbital active set, for 
2p 2 3s 4 P. NCSF is the number of CSFs in the wavefunction expansion. 



AS 


NCSF 


E 




$ sms 


ai 


a s d 


a c 


\ 








SD 


'-MR'-MCHF 








HF 


1 


-54.032 


303 


-2.5647 


3.7993 


0.7599 


3.6055 


1.5197 


3s2pld 


461 


-54.117 


540 


-2.3595 


3.5468 


0.7276 


9.0794 


1.1821 


4s3p2dlf 


2543 


-54.168 


450 


-1.3886 


3.6009 


0.7798 


5.8010 


1.1433 


5s4p3d2flg 


6447 


-54.192 


356 


-1.4432 


3.6402 


0.7761 


9.6255 


1.1132 


6s5p4d3f2g 


12210 


-54.199 


795 


-1.4511 


3.6101 


0.7584 


11.3959 


1.0968 


7s6p5d4f3g 


19832 


-54.203 


174 


-1.4498 


3.6078 


0.7589 


10.8596 


1.1146 


8s7p6d5f3g 


28014 


-54.204 


436 


-1.4505 


3.6080 


0.7615 


10.7556 


1.1106 


9s8p7d6f3g 


37595 


-54.205 


011 


-1.4509 


3.6075 


0.7615 


10.7965 


1.1068 


10s9p8d6f3g 


45958 


-54.205 


221 


-1.4505 


3.6084 


0.7607 


10.7996 


1.1094 


SDTQ-CI 


U 4s3p2dlf 


93412 


-54.205 


400 


-1.4538 


3.5997 


0.7592 


11.0883 


1.0953 


U 5s4p3d2f 


208694 


-54.205 


456 


-1.4540 


3.5977 


0.7588 


11.1538 


1.0922 


U 6s5p4d3f 


393284 


-54.205 


469 


-1.4543 


3.5972 


0.7588 


11.1705 


1.0913 



5. Comparison with experiments 

The hyperfine constants Aj and Bj are estimated for both isotopes 14 N and 15 N from the 
hyperfine structure parameters calculated with the most elaborate correlation models 
(last lines of Tables |5J E] and @J, using the nuclear data taken from Stone [T7] and 
summarized in Table |5j From equations (jSJl-dEJ), one realizes that the ratio between the 
magnetic hyperfine constants characterizing a given J-level of the two isotopes should 
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Table 5. Nuclear data for 14 N and 15 N [T7]. 



isotope 


I 


ui H 


Q (b) 


"N 


1 


+0.40376100(6) 


+0.02001(10) 


15 N 


1 

2 


-0.28318884(5) 


0.0 



be 



Aj( 15 N)/Aj( l *N) 



^( 15 N)J( 14 N) 



-1.4028. 



(15) 



/i 7 ( 14 N)/( 15 N) 

These non relativistic (NR) theoretical Aj and Bj hyperfine constants are reported 
in Table [6] and compared with those derived by Jennerich et al. [1] using saturated 
absorption spectroscopy 0. The comparison reveals huge discrepancies between theory 



Table 6. Hyperfine constant comparison between observation |pQ and theory (present 
work) : non relativistic (NR) ab initio and relativistically corrected (+RC) results. All 
values are in MHz. 



Levels 




exp. [T] 



theory 



NR 



+RC 



A 



B 



A 


B 


A 


B 


99.70 


0.0 


100.21 


0.0 


62.94 


4.10 


62.46 


4.10 


124.58 


-5.13 


124.84 


-5.12 


-53.64 


0.0 


-52.25 


0.0 


48.58 


-2.58 


51.04 


-2.95 


29.30 


3.22 


33.16 


2.57 



1/2 

P 3/2 
D 5/2 



4 no 
4 DO 

^3/2 

4 DO 

^5/2 



4 n° 

U l/2 
4 no 
^3/2 

4 n° 
A n° 

U 7/2 



+103.4(14) 
-153.1(23) a 
-47.93(48) 
-90.71(71) 

167.1(13) 
70.0(12) 
46.20(74) 

+153.1(23) 
-103.4(14) a 
92.4(17) 
41.5(14) 

-9.35(55) 



-139.85 

-88.29 
-174.75 

75.24 
-68.15 
-41.11 



-106.89 

-49.14 
-56.74 
-77.76 



-140.56 

-87.62 
-175.12 

73.29 
-71.60 
-46.52 



-104.02 

-44.49 
-51.57 
-78.04 



-69.76(90) 
+112.3(13) a 
35.52(44) 
64.76(42) 

-133.2(22) 
-48.56(74) 
-32.83(44) 

-112.3(13) 
+69.76(90) a 
-64.41(79) 
-28.19(62) 

6.31(72) 



0.0 

-0.98(48) 
-3.9(10) 

0.0 
8.69(87) 
5.0(11) 

0.0 

10.46(88) 
-0.2(15) 
-12.6(13) 



76.20 


0.0 


74.15 


0.0 


35.03 


0.0 


31.71 


0.30 


40.45 


-2.30 


36.76 


-1.69 


55.43 


-6.43 


55.63 


-6.44 



a Second proposition of Jennerich et al (see text) 



and observation. Inconsistencies appear not only in the magnitude of the parameters, 
but even in the relative signs of the parameters for the different J-levels arising from 
the same term. 



X As discussed in pQ, the hyperhne constants derived by Cangiano et al. [2\ for 4 Pj and 4 Pj agree 
qualitatively with those of Jennerich et al. PQ. 
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Rewriting equations fll])-® as 

A orb = a x Kf , (16) 

A? = l -G,g s ^a sd Kf , (17) 

Aj = l -G,g s ^ra c K c j, (18) 

Bj = -G q Qb q K'j, (19) 



the numerical factors Kf , Kj d , Kj and K'j are calculated from the expectation values 
of the angular momentum scalar products and reported in Table [7] for each J-values 
of the (odd and even) 4 P and A D terms. The consistency of the theoretical values 

Table 7. J-dependence of the hyperfine contributions to Aj and Bj constants for 4 P 
and A D terms. 



J 


R orb K sd K c Rlj 
ip 


K orb R sd K c R/j 


1/2 
3/2 
5/2 
7/2 


2 10 10 n 

3 9 27 U 
4 68 22 4 
15 45 135 5 

2 2 2 i 
5 5 15 1 


1 —I —'A 

x 3 9 u 

2 14 2 n 
5 15 45 U 

11 37 26 9 
35 105 315 25 

2 2 2 i 

7 7 21 



calculated for the different J-values within a given LS term with the relative weights of 
the different contributions making the total hyperfine constant, is obviously satisfied, by 
construction. However there is no such constraint in the experimental analysis and one 
can show that no physical set of underlying {a/, a s d, a c } parameters fits the experimental 
A j-values with the numerical factors of Table 

Some ambiguity was pointed out in the line assignment of the spectra of the 
4 Pi/2 4 -^i/2 transition for both isotopes. Due to this identification problem, 

Jennerich et al [I] determined two possible values for each of the hyperfine constants of 
the 4 Pi/2 and 4 -D°/ 2 states, both reported in Table [6j On the basis of crossover intensity 
arguments, they gave their preference to the first set. Looking to our theoretical values, 
we claim that their choice was definitely not the good one. Another surprise appears: 
they deduced a large B value for A D^ 2 14 N that contrasts with the zero numerical factor 
K' 3/2 of Table d 

6. Relativistic corrections 

Relativistic effects influence atomic wave functions in basically two ways: through 
contraction of radial orbitals and through LS term mixing. Contraction effects may be 
large for high Z, but remain relatively unimportant in light elements like Nitrogen [T8] . 
However, term mixing may affect the wave function and the computed hyperfine 
interaction constants in significant ways, especially when fine-structure levels belonging 
to different terms are closely spaced. To investigate term mixing we first perform 
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reference MCHF calculations for configuration expansions generated by SD-excitations 
from, respectively, {2s 2 2p 2 3s, 2p 4 3s} and {2s 2 2p 2 3p, 2p 4 3p} to active sets (4s3p2dlf). 
The resulting non-relativistic radial orbitals are converted to Dirac spinors using the 
Pauli approximation [19] 

P(nK; r) = P MCHF (nl; r) , (20) 

Q{nK-r) =2\-fr + -) P(riK;r) . 

Here a is the fine structure constant and k the relativistic quantum number 

K = (j + l/2)r) when l=j+rj/2, 77 = d=l . (21) 

This is followed by relativistic configuration interaction (RCI) calculations for 
configuration expansions generated by SD-excitations from the multireference sets above 
to the Dirac spinors (see ref. [201 121] for details about the conversion of radial orbitals 
and the relativistic computer codes). In relativistic theory only J is a good quantum 
number and the relativistic configuration expansions account for LS term mixing. The 
hyperfine interaction constants obtained from the MCHF and RCI wave functions are 
displayed in Table El One observes that the relativistic effects are far from negligible, 
changing the hyperfine parameters in many cases by more than 10 %. Assuming that the 
differences between the relativistic and non-relativistic values of the hyperfine constants 
in the limited calculations with active sets (4s3p2dlf) are representative of the true 
differences, we add the former to the non-relativistic results in Table [6] to obtain the 
final relativistically corrected values (in column "+RC") of the hyperfine interaction 
constants. It is clear that relativistic effects cannot explain the huge global theory- 
observation conflict. 



7. Conclusion 



The strong disagreement between theory and observation [1] is really disconcerting. 
Hyperfine parameters have indeed been estimated ab initio successfully using similar 
methods for different atomic systems such as Carbon, Nitrogen, Oxygen or Sodium, as 
described in the introduction. The present calculated hyperfine structure constants 
disagree so strongly with experiments in comparison of the achieved theoretical 
convergence of the hyperfine parameters that we presently cast doubt on Cangiano et 
al's [2] and Jennerich et a/'s pQ analysis. We therefore encourage further experimental 
spectroscopic studies and/or reinterpretation of the near-infrared spectra. Last but 
not least, the extracted isotope shifts from the same spectra revealed a significant 
unexpected J-dependence of the specific mass shifts in both multiplets. However, the 
authors themselves [I] observed that the experimental isotope shift values are critically 
dependent on the correct interpretation of the hyperfine structures of the 14 N and 15 N 
spectra. The present questioning on the experimental determination of the hyperfine 
parameters could also be relevant in their isotope shift discussion. 
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Table 8. Hyperfine interaction constants from matching RCI and MCHF calculations. 
The relativistic orbitals were obtained from the non-relativistic ones in the Pauli 
approximation. All values are in MHz. 
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85.96 
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-59.89 


0.0 
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0.0 


1.39 


0.0 
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-69.63 


-66.18 


-3.45 


49.64 - 


-2.97 


47.18 


-2.60 


2.46 


-0.37 
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^5/2 


-42.81 


-37.40 


-5.41 


30.52 


2.60 


26.66 


3.25 


3.86 


-0.65 




-102.87 


-105.74 


2.87 


73.33 


0.0 


75.38 


0.0 


-2.05 


0.0 


4 no 

U Zj2 


-40.95 


-45.60 


4.65 


29.19 


0.30 


32.51 


0.0 


-3.32 


0.30 


4 71° 

^5/2 


-48.61 


-53.78 


5.17 


34.65 - 


-1.70 


38.34 


-2.31 


-3.69 


0.61 


^7/2 


-76.34 


-76.06 


-0.28 


54.42 ■ 


-6.48 


54.22 


-6.47 


0.2 


-0.01 
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